Creation of topological states of a Bose-Einstein condensate in a plaquette 
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We study a square plaquette of four optical microtraps containing ultracold 87 Rb atoms in F = 1 
hyperfine state. In a presence of external resonant magnetic field the dipolar interactions couple 
initial mp = 1 component to other Zeeman sublevels. This process is a generalization of the 
Einstein-de Haas effect to the case when the external potential has only C4 point-symmetry. We 
observe that vortex structures appear in the initially empty uif — state. Topological properties 
of this state are determined by competition between the local axial symmetry of the individual trap 
and the discrete symmetry of the plaquette. For deep microtraps vortices are localized at individual 
sites whereas for shallow traps only one discrete vortex appears in the plaquette. States created in 
these two opposite cases have different topological properties related to C4 point-symmetry. 

PACS numbers: 03.75.Mn, 05.30.Jp, 75.45.+j 



I. INTRODUCTION 

Experimental achievement in cooling alkali atoms pro- 
vide a multiple possibilities to explore Bose-Einstein con- 
densation phenomena. The first Bose-Einstein conden- 
sates were created in magnetic traps and allowed for 
studying atomic condensates properties related to weak 
contact interactions. Soon after that an optical confine- 
ment of ultracold atoms became possible and new sys- 
tems - spinor condensates - appeared. Here, the ability 
to simultaneously trap atoms in many magnetic states 
makes it possible to investigate nontrivial spin dynamics 
and magnetism in ultracold Bose gases already studied 
both experimentally and theoretically QUI]. 

It was realized quite early that dipolar interactions be- 
tween ultracold species might lead to many novel phe- 
nomena • The interesting features of dipolar interac- 
tions are rooted in their long-range as well as anisotropic 
character. The Bose-Einstein condensate in 52 Cr atoms, 
whose magnetic dipole moment equals 6fj,B was obtained 
after many years of experimental struggles [l(| E| • This 
achievement opened a possibility to stud y a variety of 
phenomena related to dipolar interactions |l2j.[l3|. Some 
authors noticed however, that magnetic dipole 

interactions might lead to observable effects even in a gas 
of rubidium atoms. An example is the Einstein-de Haas 
effect [13, EH (studied also for 52 Cr condensate in [l9j). 
It's main idea is that rubidium atoms trapped initially 
in mp = +1 component can be efficiently transferred to 
mp = state and then further to mp = —1 Zeeman sub- 
level. This can happen only if external magnetic field is 
tuned to a resonant value. Moreover, atoms in mp = 
and mp = — 1 states acquire one or two quanta of orbital 
angular momentum, respectively (T7| . 

Ability of making use of interfering laser beams to 
create optical lattices filled with ultracold bosonic or 
fermionic atoms opened new possibilities. Atomic physics 
started to penetrate areas traditionally associated with 
solid state physics. Observation of the Mott insulator 
- superfluid transition in an optical lattice [20M22I ] was 



a very important step towards investigation of strongly 
correlated systems described by various Hubbard-type 
Hamiltonians. Since then degenerate gases inside opti- 
cal lattices are studied very intensively [23L l24j . 

The progress in experimental techniques allowed to 
create new types of optical lattices such as superlattices 
or arrays of coupled plaquettes [25[ . In this article we fo- 
cus on topological states of the ultracold atomic system 
in a single plaquette. We study a case of relatively large 
occupation of the plaquette what allows for the mean 
field description of the system. Such systems are ex- 
perimentally accessible [2^, [27j • We concentrate on the 
topological phases of the macroscopic order parameter. 
Creation of these dynamical states implicitly utilizes the 
dipolar interactions tuned by the external magnetic field. 
External magnetic field is a perfect tool for manipulat- 
ing the system in the controlled manner. The analysis 
presented in this article becomes important especially in 
the context of recently performed experiment (27j aim- 
ing at the observation of the Einstein-de Haas effect in a 
chromium condensate in an optical lattice. 

If the atoms are located in an array of traps then the 
Einstein-de Haas effect may exhibit a new interesting as- 
pects. Conservation of the total angular momentum is 
the essence of the Einstein-de Haas effect. In this article 
we assume that 87 Rb atoms are placed in a square plaque- 
tte of four optical traps. Therefore, the external potential 
has no axial symmetry and z-component of the total an- 
gular momentum has not to be conserved. Instead, the 
C4 point-symmetry comes into a play. The aim of this 
article is to investigate the role of this symmetry in the 
dynamics driven by the dipolar interactions. The relative 
weight of the local axial symmetry of the trap and the 
discrete symmetry of the plaquette can be controlled by 
varying the height of the barrier separating the traps. 

This paper is organized as follows. In Sec. |H] we de- 
scribe the system under consideration. In particular, we 
discuss in detail the properties of dipolar interactions 
which are essential for the physics reported in the ar- 
ticle. Sec. IIHI presents numerical results regarding the 
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Einstein-de Haas effect for a rubidium condensate in a 
plaquette whereas Sec. IIVI explains why dipolar reso- 
nances sometimes lead to the global discrete vortex state 
while otherwise to the array of local vortices. In Sec. |V| 
we discuss the stability of states obtained via the reso- 
nant Einstein-de Haas effect. We end with conclusions in 
Sec. EH 



position r are defined as: 



II. DESCRIPTION OF THE SYSTEM 



The single-particle Hamiltonian has the following form: 
H = | dr ft(r) + V tr (r) - 7 Bf) $(r) , (1) 

where -0(r) = (ifj + (r), ipo( r ), '4 ; -( r )) T is the spinor anni- 
hilation operator. F = (F x , F y , F z ) are standard F = 1 
spin matrices, and 7 is the gyromagnetic coefficient. The 
external potential of the plaquette Vt r is: 



V tr (r) = V (cos 2 {k x) + cos 2 (fc y)) 
+ ^^Kt^ + V 2 ) + u z z 2 ) . 



(2) 



The first term describes periodic lattice while the second 
one corresponds to a relatively weak harmonic potential. 
The harmonic confinement is essential when the stability 
of topological states created by resonant Einstein-de Haas 
effect is discussed (Sec. IV]) . Otherwise (Sees. Mil and 
UV|) . we assume that only axial (i.e., along z-axis) part 
of harmonic trapping is present. The lattice is defined 
on a square of the xy plane centered at x = y = and 
containing four lattice minima which form a square and 
are separated by c? = tt/^o- Vo is a potential barrier 
separating lattice sites. The last term in (JXJ) is responsible 
for the linear Zeeman shift due to the uniform magnetic 
field B directed along the z-axis. In the following we 
assume that the magnetic field is weak and therefore the 
quadratic Zeeman effect can be neglected. 

The short-range interactions are typically described by 
a pseudo-potential. In the case of spinor F = 1 con- 
densate this interaction can be splited into the spin- 
independent term proportional to Co = 4irti 2 (2a,2 + 
ao)/(3m) and the spin-dependent part proportional to 
C2 = 47rfi, 2 (a2 — ao)/(3m) [28( where 02 (ao) is a scatter- 
ing length of two colliding atoms with total spin F = 2 
(F = 0). The Hamiltonian describing the short-range 
interactions can be written in the form: 



Hr 



dr 



n 2 (r) 



£2 
2 



(3) 



where n(r) = V'lV's is the total atom density, :: de- 
notes normal ordering, and square of the total spin op- 
erator can be written as: F 2 (r) = 2 (F + (r)F_(r) + 
F_(r)F + (r)) + F z (r)F z (r). The raising (lowering) op- 
erators of the z-component of the spin of the atom at 



F ± (r) = ^t (r) ^k±^k^ (r)) 



(4) 



while the magnetization density is: 

F z (T)=ft(T)F g #(T) . (5) 

Finally, the long-range dipolar Hamiltonian is: 



H n = %r I dr / dr 



, _ F(r)F(r') - 3 [F(r)n] [F(r')n] 



| r _ r /|3 

(6) 

where n is a unit vector in the direction of r — r'. Using 
explicit form of spin-1 matrices the above Hamiltonian 
can be brought to the form: 

ffD = ^/ dr / dr '^F : ^ (r ' r ' ): ' (7) 

where /i£>(r,r') has the form: 

h D = A (F z (r')F z (r) - F + (r')F_(r) — F_(r')F + (r)) 

- 3 sin 2 6 '{e- 2 ^F + {v')F + {T) + e 2i *F_ (r')-P_ (r)) 
- 3/2sin20e-^( J F + (r')F 2 (r)+ J F 1 z (r')F + (r)) 

- 3/2sin26»e I *(F 2 (r') J F-(r)+F_(r')i ;l 2 (r)). (8) 

9 and (j> are the spherical angles of the vector r — r' con- 
necting two interacting atoms and A = 1 — 3 cos 2 0. The 
form of ([5]) allows for physical interpretation of all terms. 
The first line represents dipolar interactions which do 
not lead to change of total magnetization of the field: 
^-components of spin of both interacting atoms remain 
unchanged or the z-component of one atom decreases by 
one while the z-component of the second atom increases 
by one. The second line collects terms describing pro- 
cesses where both interacting atoms simultaneously flip 
the z-axis projection of the spin: both by +1 or both by 
— 1. Notice that corresponding terms are multiplied by 
the phase factor e -2 ^ or e 2i< *, respectively. This corre- 
sponds to change of the projection of the orbital angular 
momentum of atoms in their center of mass frame by 
—2 or 2 quanta. The last two lines describe processes 
in which the spin of one interacting atom is unchanged 
while the z-axis component of the spin of the other atom 
changes by +1 (or by —1). This spin flipping term is 
multiplied by the phase factor e"*^ (or e^) what signifies 
the change of the z-projection of relative orbital angular 
momentum of interacting atoms by —1 (or +1 respec- 
tively). Evidently, the dipolar interactions conserve the 
z-projection of the total angular momentum of interact- 
ing atoms. 



III. EINSTEIN-DE HAAS EFFECT IN A 
PLAQUETTE 



We study a system of N 
tially in the polarized mp = 



4 x 10 3 87 Rb atoms ini- 
+1 state and the ground 
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state of the external potential. Because occupation of 
this state is macroscopic we describe the system by 
a spinor wavefunction satisfying the Gross-Pitaevskii 
(GP) equation. We therefore replace the field opera- 
tor ip(r) by the c- number spinor wavefunction ip(r) = 
(■0-1- (r), ipo( r ), i 1 - ( r )) T - Standard optical lattices for ru- 
bidium atoms have the site separation in the submicron 
regime. In such a geometry and the number of atoms con- 
sidered here the tunneling is large even for barriers much 
larger than it is accessible experimentally. However, it 
is possible to suppress the tunneling by using larger sep- 
arations between microtraps. Large separations can be 
reached in various ways, for example by varying the angle 
between co-propagating laser beams creating a periodic 
potential. 

We solve the spinor GP equation with the nonlocal 
dipolar potential on a spatial grid 32 x 32 x 16 points. 
The spatial steps are dx = dy = 0.25/im and dz = 0.5^m, 
the wavevector is ko = 0.74//im, and the axial confine- 
ment is given by Hlj z /Eji = 0.32, where En = fi, 2 fcg/2rn 
is the single-photon recoil energy. Note that at each site, 
the microtrap potential has a prolate shape due to the 
tight confinement in the xy plane in comparison with the 
one along the z-direction according of the harmonic trap. 
Here, we assume that the radial harmonic confinement is 
absent (luj_ = 0). We instantly turn on the magnetic 
field and monitor dynamics of the spinor wavefunction. 
Transfer of atoms to other Zeeman states strongly de- 
pends on the value of the external magnetic field and 
geometry of the trap [l7|. This transfer is possible only 
when the Zeeman energy compensates for the excitation 
energy. Therefore we have to adjust the value of the ex- 
ternal magnetic field to the resonant value. For a given 
initial state there are many magnetic resonances (29| . We 
focus on the first resonance, corresponding to the smallest 
value of the magnetic field and the most efficient trans- 
fer of atoms to other Zeeman components. We start our 
simulations with relatively low barrier between plaquette 
sites, Vo = lOEfj. After about t = 4s, at resonant mag- 
netic field B = 0.075 mG, a transfer of atoms to mp = 
component reaches maximum and starts do decrease. We 
observe a kind of oscillations of population at longer time 
scale if magnetic field is on. To check the stability of the 
created dynamically state we switch off the magnetic field 
at the moment when the transfer is maximal. This way 
the dipolar interactions are tuned out of the resonance. 
Also the spin-mixing contact term is small because the 
number of atoms in mj? =0 component is of the order 
of 100 only. Therefore further dynamics looks stationary 
at least at a time scales of few seconds. Population of 
mp = —1 component remains always very small, how- 
ever its topological structure proves dipolar character of 
atoms transfer. Analysis of a one-particle density matrix 
averaged over z-axis (30| indicates that all three compo- 
nents of the spinor wavefunction are coherent. In the 
following discussion we ignore mp = — 1 component be- 
cause of its negligible occupation. It is, however, present 
in our calculations. 
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FIG. 1: Phase (left), the density (middle) in the xy plane at 
z = 1.3/fco and the density in the xz plane at y — 1.5/fco 
(right) of ttif — component at t = 4 s for Vo = IOEr in 
the case without radial harmonic confinement. Here, fco = 
0.74/^m and the distance between lattice sites equals do = 
Ti/ko = 4.24/xm. Total number of atoms N = 4 x 10 3 and the 
resonant magnetic field B — 0.075 mG. Maximal transfer to 
niF = component is equal to No = 50 atoms at t — 4 s. 



In Fig. [T] we present the density distribution and the 
phase oi mp = component. Atomic cloud is con- 
centrated around two horizontal xy planes, the one at 
z = 1.3/fco & n d the second at z = —1.3/fco- I n both 
planes the density is the same. Discrete symmetry of the 
plaquette is clearly visible - almost all atoms are located 
at centers of the quadrants of the plaquette. The phase 
of the wavefunction is constant in each quadrant but is 
not constant over the entire plaquette: it jumps by n/2 
between neighboring quadrants. As a result the phase 
of the atomic wavefunction winds up by 27r around the 
plaquette center. Evidently, the singly quantized discrete 
vortex in mp — component is present. The geometry of 
the final state proves that the discrete symmetry of the 
whole plaquette determines the dynamics of the Einstein- 
de Haas effect. It dominates the local axial symmetry of 
the individual site in this case. 
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FIG. 2: Transfer of atoms to mp = state as a function 
of magnetic field for the plaquette with Vo = WEr. Two 
resonances are visible. The one for lower magnetic field cor- 
responds to the case when the transfer of atoms is governed 
by the discrete symmetry of the plaquette (see Fig. [TJ. For 
the second resonance (for higher magnetic field) the transfer 
to my =0 state occurs locally, at individual plaquette sites 
(see Fig. S}. 
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However, the topology of the final state might change. 
It turns out that there exists a resonance which populates 
the state with local vortices located at the lattice sites. 
To find such a resonance we still keep the height of the 
barrier separating the sites at Vq = IOEji and increase 
the magnetic field. We observe a significant transfer of 
atoms to mp = state for the magnetic field equal to 
B = 0.26 mG (see Fig. [2]). In Fig. [3] we show the density 
distribution of mp = component at the maximum of 
transfer. Similarly as in the previous case the final state 
does not change noticeably when evolved with external 
magnetic field set to zero. The number of mp = atoms 
in each site is equal to No = 180. We see the array of 
singly quantized vortices. At each lattice sites the den- 
sity forms two rings - the one above z = plane and the 
second below this plane. The phase of the corresponding 
wavefunction winds up by 2tt around every site center. 
Even more, the phases of individual vortices seem to be 
correlated because lines where the phase is zero are par- 
allel. Evidently, in the case of this resonance the transfer 
of atoms to mp = component takes place locally, at 
every lattice site and resembles the Einstein-de Haas ef- 
fect in axially symmetric harmonic potential. Note that 
the potential in our case is axially symmetric only close 
to the minimum at every site of the plaquette. 




FIG. 3: Phase (left), the density (middle) in the xy plane at 
z — 1.3/fco and the density in the xz plane at y = 1.3/fco 
(right) of m_F = component at t = 4 s for Vo = 10 En in 
the case without radial harmonic confinement. Here, fco = 
0.74/^m and the distance between lattice sites equals do = 
Ti/ko = 4.24/xm. Total number of atoms N = 4 x 10 3 and the 
resonant magnetic field B = 0.26 mG. Maximal transfer to 
771 f = component is equal to No = 180 atoms at t = 0.6 s. 



IV. APPEARANCE OF DISCRETE VORTEX 
AND LOCAL VORTICES STATES 

To get better understanding of our numerical results 
we will limit our spinor wavefunction to a few relevant 
states only. At every site localized at = (xi , j/j, 0) 
where i — 1, . . . , 4 (for 2x2 plaquette) we keep in the 
mp = 1 component only the ground state: 

(r\i) g =N' Wo(x-x i ,y-y i )e-^^. (9) 

Two-body dipolar interactions can transfer only one par- 
ticle from the initial ground state (r\i) g to excited states 
with mp = 0. Simultaneous transfer of both particles is 



strongly suppressed for large initial state occupation (2j| . 
There are two relevant low energy excited states (r|i) ei 
and (r|i) e2 in a given site which can be resonantly pop- 
ulated. The first one is the state with no orbital angular 
momentum and one z-excitation quantum: 

(r|i) ei =Af 1 W {x-x i ,y-y i )ze-^/ 4 , (10) 

while the second one has one quantum of orbital angular 
momentum and one z-excitation quantum: 

(r|i) ea =Af 2 W 1 (x-x i ,y-y i )ze- 0z2 / A . (11) 

Wo(x,y) is the 2D ground Wannier wavefunction with 
zero orbital angular momentum being a product of ID 
ground Wannier states Wo(x,y) = Wo(x)Wo(y) while 
Wi(x,y) = [W 1 (x)Wo(y)+iW (x)W 1 (y)]/V2is the 2D 
excited Wannier state with one quantum of orbital an- 
gular momentum along the z-axis. Here W\ is the first 
excited ID Wannier state, j3 = hu> z /En is the harmonic 
oscillator energy expressed in units of the recoil energy 
while Ao, A/i, N% are normalization constants. For large 
occupation of the initial state the processes leading to 
the chosen excited states strongly dominate others as dis- 
cussed in [2{j. We will neglect the mp = — 1 component 
because its occupation is due to the second order per- 
turbation in dipolar interactions. The first excited state 
(JTUJ) could, in principle, be responsible for appearance of 
discrete vortex in the entire plaquette, while the second 
one (fTTjl corresponds to the vortex array. 

Mean energies of all above states do not depend on the 
lattice site and they are E g , E ei and E e2 respectively. 
All states, i.e. the ground and the two excited states are 
four-fold degenerate in the 2x2 plaquette. Tunneling 
however couples neighboring sites and partially removes 
this degeneracy. In the chosen basis the spectrum of the 
single particle Hamiltonian splits into a three 'bands' of 
four Bloch states: 



\i)a = 




+ |2) a - 


H3)a- 


H4>a), 


(12a) 


\0)a = 




+ t\2) a 


"|3>a 


-i|4)„), 


(12b) 


\0)a = 




-i\2) a 


"|3>a 


+ *|4)a), 


(12c) 


lt)« = 


|(H>- 


"|2)a- 


H3)a- 


-|4>a), 


(12d) 



where a = g, ei, e-2 an d the state vector \i) a corresponds 
to the wavefunction (r\i) a . Energies of the four states 
within a given band are equal to {E a — 2 J a , E a , E a , E a + 
2J a }, respectively. Tunneling coefficients in the ground 
band and in the first excited bound are positive and equal 
J g = J ei because all states in these bands have the same 
spatial profile in the x-y plane. Therefore the states \\.) g 
and |4) ei are the lowest energy states in corresponding 
bands. 

Tunneling between vortex states in the second excited 
band is negative J e2 < 0. Therefore the lowest energy 
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state in this band is anti- ferromagnetic |t) e2 with alter- 
nating relative phases rather then the ferromagnetic one 
1 4-}e 2 - Subspace of eigenenergy E a spanned by vectors 
\0) a and \0) a is two-dimensional, so one can choose any 
combination of these vectors as the basis. We choose 
these particular combination since such vectors are not 
only the eigenstates of the Hamiltonian but also they are 
the eigenstates belonging to the proper representations of 
the C4 symmetry group with nonzero eigenvalues ±1, i.e., 
they are closely related to the eigenstates of the angular 
momentum operator in the continuous case. 

In our numerical simulations we prepared the system in 
the \\.} g state. Dipolar interactions can couple this state 
to the both excited bands. The strength of this coupling 
is given by the two-body dipolar matrix element between 
the initial and final two particle states. In the final state 
one atom remains in the initial state but the second one 
occupies one of the states from the excited bands. Direct 
evaluation of the relevant integrals shows that only two of 
them strongly (by several orders of magnitude) dominate 
the others. They correspond to the process when one 
atom is transferred to the discrete vortex state in the 
first excited band, i.e., the state in which density profile 
is almost the same as in the initial state but the phase 
changes from site to site by tt/2: 

\i) g -> I0) ei = \ (|l) ei + i|2) ei - |3) B1 - i|4) ei ) , (13) 
or to the state with the array of four vortices: 

\i)g -> \i)e 2 = \ (|l)e 2 + |2) e2 + |3) e2 + |4) e J . (14) 

The energy difference for the first process is SEi = E ei — 
Eg + 2J g while for the second one is SE 2 = (E e2 — E g ) — 
2(J e2 — Jg). Evidently resonant value of the magnetic 
field for the first process is smaller then for the second 
one. Moreover, in the case of the second resonance the 
relative phases of vortices created at individual lattice 
sites are the same as relative phases in the initial state. 
The phase between vortices at different lattice sites is 
inherited from the initial state. If the initial state is \\) g 
then atoms are transferred to the analogous state in the 
second excited band, |t)e 2 - This observation agrees very 
well with our numerical results presented in the previous 
section. 

Now it is clear that due to the dipolar interactions 
transfer of angular momentum from the spin part to the 
spatial part can be realized in two ways: by exciting par- 
ticles in every site to the second excited band with array 
of vortices (e?) or by transferring particles to the first ex- 
cited band - not changing the local density but changing 
a relative phases between plaquette sites (ei). 

Exact values of these two relevant dipolar transition 
amplitudes obviously depend on the barrier height and 
the strength of confinement in z-direction (parameter /3). 
Since for deeply enough lattice the sites become more 
symmetric, transition to on-site vortices should dominate 
over the transition to discrete vortex in the plaquette. 
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FIG. 4: Transition amplitudes induced by dipolar interactions 
for the 2x2 plaquette with hoj z /En = 3. For deeper lattices 
transition to the four independent vortices located at lattice 
sites (solid line) always dominates over the transition to the 
discrete vortex (dotted line). 

From the other side, when the lattice is shallow the on- 
site Wannier states with vortex are quite different from 
eigenfunctions of 'local' (at a given site) orbital angular 
momentum operators. Moreover the state with discrete 
vortex, due to a high tunneling, becomes quite almost in- 
variant under rotation around the center of the plaquette 
by arbitrary angle (up to the phase factor) . Therefore we 
expect that in this case transition to the plaquette vor- 
tex will be dominant. We have confirmed all these pre- 
dictions by numerical calculation of dipolar amplitudes. 
The results are presented in Fig. 0J The height of the 
barrier for which two considered amplitudes are equal 
as well as transition amplitudes depend on the parame- 
ter (3 characterizing the trapping potential. We should 
emphasize again that amplitudes which initially do not 
conserve total angular momentum, i.e., the ones for final 
states with opposite or without vortex, are numerically 
many orders smaller than considered ones. 

Presented model is in a very good qualitative agree- 
ment with the results of numerical experiment presented 
here. However, the values of magnetic field when pop- 
ulation of the discrete vortex dominates population of 
vortex array state is different. This is because of simpli- 
fications we made to get a clear picture of relevant pro- 
cesses. In particular the contact interactions were totally 
ignored. Their role is to change both energies and spatial 
spreading of the atomic states what inevitably leads to 
quantitative differences. 

V. STABILITY OF TOPOLOGICAL STATES 

Now, we are going to investigate the stability of states 
populated by resonant Einstein-de Haas effect. For that, 
we must turn on the radial harmonic confinement since 
studying the stability of discrete vortex and local vortices 
states we will first remove the optical potential and next 
restore it. 

We first checked that in the case when the radial har- 
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monic confinement is present together with the periodic 
potential, we are able to find the same kinds of reso- 
nances as m Sec. [Till Indeed, the presence of an extra 
radial harmonic potential results in a small shift of the 
resonant magnetic field and a small change in a num- 
ber of atoms transformed to the mp = state as com- 
pared to the cases discussed in Sec. IIIII For example, 
for the optical trap with Vq = IOEr and the resonance 
which populates the singly quantized discrete vortex in 
mp = component we find that the resonance magnetic 
field equals B = 0.08 mG and the number of transferred 
atoms is A*o = 80. For further studies we focus on the 
resonance just mentioned and the other one which leads 
to significant occupation of vortex states at individual 
plaquette sites and is found for parameters Vq = 2QEr, 



0.85 mG, and N = 160. 
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FIG. 5: Sequence of images of the phase (upper panel) and 
the density (lower panel) in xy plane (z = 1.3/fco) of tuf = 
component during lowering and raising of the barrier sepa- 
rating microtraps: Vo = WEr (left), Vq = OEr (middle), 
and Vo = WEr (right). The barrier was lowered linearly in 
time during t = 0.5s. The images in the middle column are 
shown at t = 0.15s after the barrier is removed. Then the 
barrier was raised in t = Is and the initial configuration was 
recovered (see right column). 




FIG. 6: Sequence of images as in Fig. [S] but for Vo = 
20Er,0Er and 2QEr (left, middle, and right column, re- 
spectively). The images in the middle column are shown at 
t = 0.15s after the barrier is removed. Then the barrier was 
raised in t = Is and the initial configuration was recovered 
(see right column). 



relatively low density surrounding the first one can be 
seen, Fig. [HI However, after raising back the barrier both 
initial configurations are recovered (see the right columns 
of Figs.[5j[6|) what indicates that the initial states are es- 
sentially different and one cannot switch between these 
states after they are created. 
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Depending on the barrier separating lattice sites and 
the value of an external magnetic field the resonant dipo- 
lar interactions drive the system towards the states of 
a different geometry. It is interesting to check if these 
states are robust against perturbation of the trapping 
potential of the lattice. To this end we ramp adiabati- 
cally the lattice potential down to zero and then we rise 
it back to the initial value for both studied cases. The 
sequence of pictures in Figs. [5l [6] shows how the density 
and the phase of tof = component change while the 
microtrap barrier is first lowered and then raised back. 
The lowering process lasts for t = 0.5 s. The middle 
columns show the density and the phase at 0.15 s after 
the barrier height reached zero. We see that in both cases 
atoms form a vortex around the center of the plaquette. 
In the case of high initial barrier an additional ring of 



FIG. 7: Phase portraits of ttif =0 component obtained at 
Vo = 20Er (left panel) and Vo = WEr (right panel). Notice 
a singly quantized discrete vortex of positive charge corre- 
sponding to Vo = WEr case. For Vo = 20Er in addition to 
four singly quantized vortices of positive topological charge lo- 
cated at the center of each quadrant where the atomic density 
is large, there are five more vortices. One of them is located 
at the plaquette center and the other four have opposite topo- 
logical charge and are located at the bonds connecting lattice 
minima. These additional five vortices are hidden in the re- 
gion of a very low atomic density. The vortices and their 
charges are marked by white and black circles. 



To understand this difference let us look once more into 
the phases of both states as they were created via dipo- 
lar interactions, Fig. [7J As the discrete vortex, obtained 
for lower barrier, Vq = 10_Er, is the only singularity of 
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the phase of the mp — component (right panel), much 
richer structure can be observed for the array of vortices 
obtained at Vb = 20-Er (left panel). In addition to four 
singly quantized vortices of positive topological charge 
located at the center of each minimum there is one more 
vortex of the same charge located at the plaquette cen- 
ter and four more vortices of opposite topological charge. 
These vortices and their charges are marked by white and 
black circles in Fig. [7J The additional five vortices are 
hidden in the region of a very low atomic density. The 
number of vortices strongly differentiates between these 
states. The vortex structures are stable against perturba- 
tion of the lattice potential therefore it is not possible to 
switch between them by changing the intertrap barrier. 
In Fig. [6] only one vortex is clearly seen. Other vortices 
are located at regions of vanishing density. However, 
numerical evaluation of circulation around every point 
of the grid allows to detect them. The system recovers 
the initial symmetry even after the barrier was removed. 
This memory effect is due to different number of vortices 
in both studied states. These two states belong to par- 
ticular representation of the C4 point-symmetry group. 
This representation is characterized by the presence of 
singly quantized vortex at the center of the plaquette. It 
is worth to note at this point that an extensive discussion 
of solutions of the GP equation belonging to different rep- 
resentations of C n point-symmetry group can be found 
in Ref. HH]. 

We would like to stress that our findings might be im- 
portant for the observation of the Einstein-de Haas effect 
in optical lattices. Some experimental effort has been 
already taken in this direction [27j. Our results show 
that the properties of the rotating states created via the 
Einstein-de Haas effect strongly depend on the geometry 



of the trapping potential and the value of the external 
magnetic field. 



VI. CONCLUSIONS 

In conclusion, we studied formation of topological 
states via resonant dipolar interactions. We paid a par- 
ticular attention to competition between the local axial 
symmetry of the individual microtrap and global discrete 
symmetry of the entire plaquette. We showed that the 
role of these two symmetries can vary depending on the 
height of the intertrap barrier. For shallow enough pla- 
quettes when the mean field wavefunction penetrates the 
whole plaquette a single discrete vortex is created. In the 
opposite case of low tunneling the transfer of atoms has a 
local character and axial symmetry of individual site al- 
lows for formation of the array of vortices located at each 
plaquette site. However, for a wide range of intermediate 
lattice's depths both kinds of topological states can be 
populated depending on the value of applied magnetic 
field. These topological states are created dynamically 
with the help of resonant Einstein-de Haas effect. More- 
over, the vortex structures present in both states make 
these states very stable against perturbation of the lat- 
tice potential. Even when the discrete symmetry of the 
plaquette is destroyed by removing the lattice potential, 
the system comes back to the initial configuration after 
the lattice potential is raised back. 
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